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^ ■ ABSTRACT 

1) , The strong coupling as{M'^) has been measured using hadronic decays of bosons 

Qj ' collected by the SLD experiment at SLAC. The data were compared with QCD predic- 



tions both at fixed order, O^aj), and including resummed analytic formulae based on 



^ . the next-to-leading logarithm approximation. In this comprehensive analysis we stud- 

H ; 

5t , ied event shapes, jet rates, particle correlations, and angular energy flow, and checked 

the consistency between q;s(M|) values extracted from these different measures. Com- 
bining all resuhs we obtain a^(M|) = 0.1200 ± 0.0025(exp.) ± 0.0078(theor.), where 
the dominant uncertainty is from uncalculated higher order contributions. 
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1. Introduction 



Achieving precision tests of the Standard Model of elementary particle interactions 
is one of the key aims of experimental high energy physics experiments. Some mea- 
surements in the electroweak sector have reached a precision of better than 1% 
However, measurements of strong interactions, and hence tests of the theory of Quan- 
tum Chromodynamics (QCD) p|, have not yet achieved the same level of precision. 
This is largely due to the difficulty of performing QCD calculations, both at high order 
in perturbation theory and in the non-perturbative regime, where effects due to the 
hadronization process are important. QCD is a theory with only one free parameter, 
the strong coupling a^, which can be written in terms of a scale parameter Aj-fg. All 
tests of QCD can therefore be reduced to a comparison of measurements of a^, either 
in different hard processes, such as hadron-hadron collisions or e"'"e~ annihilations, or 
at different energy scales Q. In this paper we present measurements of as in hadronic 
decays of bosons produced by e^e~ annihilations at the SLAC Linear Collider and 
recorded in the SLC Large Detector. 

Complications arise in making accurate QCD predictions. In practice, because of 
the large number of Feynman diagrams involved, QCD calculations are only possible 
with present techniques to low order in perturbation theory. Perturbative calculations 
are performed within a particular renormalization scheme 0, which also defines the 
strong coupling. Translation between different schemes is possible, without changing 
the final predictions, by appropriate redefinition of as and of the renormalization scale 
1^. This leads to a scheme- dependence of a^, which can be alleviated in practice by 
choosing one particular scheme as a standard and translating all as measurements to 
it. The modified minimal subtraction scheme (MS scheme) is presently used widely 
as this standard. An additional complication is the truncation of the perturbative 
series at finite order, which yields a residual dependence on the renormalization scale, 
often denoted by /i or equivalently by / = /x^/Q^, which then becomes an arbitrary 
unphysical parameter. 
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In our previous studies of jet rates and energy-energy correlations ^ it was 
shown that the dominant uncertainty in ^^(Mf ) measurements arises from this renor- 
malization scale ambiguity. Given that infinite order perturbative QCD calculations 
would be independent of /i, the scale uncertainty inherent in measurements is a 
reflection of the neglected higher order terms. 

Distributions of observables in the process e^e~ hadrons have been calculated 
exactly up to 0{a1) in QCD perturbation theory 0. One expects a priori that the 
size of the uncalculated 0{al) and higher order terms will in general be different for 
each observable, and hence that the scale dependence of the values measured using 
different observables will also be different. In order to make a realistic determination 
of and its associated theoretical uncertainty using 0{a1) calculations it is therefore 
advantageous to employ as many different observables as possible. Our previous mea- 
surements of as(M|) were based on extensive studies of jet rates and energy-energy 
correlations and their asymmetry using approximately 10,000 hadronic decays 
collected by the SLD experiment in 1992. In this comprehensive analysis we have used 
the combined 1992 and 1993 data samples, comprising approximately 60,000 events, 
to make an improved determination of as{M^) using fifteen observables presently cal- 
culated up to O^a^) in perturbative QCD. 

In addition, for six of these fifteen observables, improved calculations can be for- 
mulated incorporating the resummation []TU|, [TT], |T^, |T^, |T3[ of leading and next- 
to-leading logarithms matched to the 0{a^) results; these matched calculations are 
expected a priori both to describe the data in a larger region of phase space than the 
fixed-order results, and to yield a reduced dependence of on the renormalization 
scale. We have employed the matched calculations for all six observables to determine 
q;s(M|), and have studied the uncertainties involved in the matching procedure. We 
have compared our results with our previous measurements and with similar measure- 
ments from LEP. 

We describe the detector and the event trigger and selection criteria applied to the 
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data in Section 2. In Section 3 we define the observables used to determine q;s(M|) in 
this analysis. The QCD predictions are discussed in Section 4. The analysis of the data 
is described in Section 5, and a summary and conclusions are presented in Section 6. 



2. Apparatus and Hadronic Event Selection 

The e^e~ annihilation events produced at the resonance by the SLAG Linear Col- 
lider (SLC) have been recorded using the SLC Large Detector (SLD). A general de- 



scription of the SLD can be found elsewhere [|T6[. Charged tracks are measured in 
the central drift chamber (CDC) and in the vertex detector (VXD) |T^. Momentum 
measurement is provided by a uniform axial magnetic field of 0.6 T. Particle ener- 
gies are measured in the Liquid Argon Calorimeter (LAC) [|1^, which contains both 
electromagnetic and hadronic sections, and in the Warm Iron Calorimeter ]TP|. 



Three triggers were used for hadronic events. In the 1993 (1992) runs the first 
required a total LAC electromagnetic energy greater than 12 GeV (8 GeV); the second 
required at least two well-separated tracks in the CDC; and the third required at least 
4 GeV (8 GeV) in the LAC and one track in the CDC. A selection of hadronic events 
was then made by two independent methods, one based on the topology of energy 
depositions in the calorimeters, the other on the number and topology of charged 
tracks measured in the CDC. 

The analysis presented here used the charged tracks measured in the CDC and 
VXD. A set of cuts was applied to the data to select well-measured tracks and events 
well-contained within the detector acceptance. The charged tracks were required to 
have (i) a closest approach transverse to the beam axis within 5 cm, and within 10 cm 
along the axis from the measured interaction point; (ii) a polar angle Q with respect to 
the beam axis within | cos^ |< 0.80; and (iii) a momentum transverse to the beam axis, 
p_L > 0.15 GeV/c. Events were required to have (i) a minimum of five such tracks; (ii) 
a thrust axis ||2^ direction within | cos 6t \< 0.71; and (iii) a total visible energy Ems 
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of at least 20 GeV, which was calculated from the selected tracks assigned the charged 
pion mass. From our 1992 and 1993 data samples 37,226 events passed these cuts. 
The efficiency for selecting hadronic events satisfying the | cos 6t \ cut was estimated 
to be above 96%. The background in the selected event sample was estimated to be 
0.3 ±0.1%, dominated by t^t~ events. Distributions of single particle and event 

topology observables in the selected events were found to be well described by Monte 
Carlo models of hadronic decays [pl|, combined with a simulation of the SLD. 



3. Definition of the Observables 

In this section we present the definitions of the quantities used in our measurement of 
q;s(M|). We used observables for which complete 0{a1) perturbative QCD calculations 
exist. These include six event shapes, jet rates defined by six schemes, two particle 
correlations, and an angular energy fiow. 

A. Event Shapes 

Various inclusive observables have been proposed to describe the shapes of hadronic 
events in e"'"e~ annihilations. We considered those observables which are coUinear and 
infrared-safe, and which can hence be calculated in perturbative QCD. 
Thrust T is defined by Q 



T = max 





P 






Pi 





where pi is the momentum vector of particle i, and ut is the thrust axis to be de- 
termined. We define r = 1 — T. For back-to-back two-parton final states r is zero, 
while < r < I for planar three-parton final states. Spherical events have t = |. An 
axis TijYiQ^j can be found to maximize the momentum sum transverse to ut- Finally, an 
axis fimin is defined to be perpendicular to the two axes n-r and n^aj- The variables 
thrust-major T^aj and thrust-minor T^m are obtained by replacing in Eq. (|l]) by 
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nmaj or fimin, respectively. The oblateness O is then defined by 

O T^dj T^in- (2) 

The value of O is zero for coUinear or cylindrically symmetric final states, and extends 
from zero to 1/ -\/3 for three-parton final states. 

The C-parameter is derived from the eigenvalues of the infrared-safe momentum 



tensor ||24|| : 





Pi 1 




Vi 





e,. = ^^::"!^jr , (3) 



where is the p-th component of the three momentum of particle and i runs over 
all the final state particles. The tensor 6'pcr is normalized to have unit trace, and the 
C-parameter is defined by: 

C = 3(AiA2 + A2A3 + A3Ai), (4) 

where Aj (z = 1,2,3) are the eigenvalues of the tensor Opij. For back-to-back two- 
parton final states C is zero, while for planar three-parton final states < C < 2/3. 
For spherical events C = 1. 

Events can be divided into two hemispheres, a and b, by a plane perpendicular to 
the thrust axis ut- The heavy jet mass Mh is then defined as 



Mh = max{Ma,Mb), (5) 

where Ma and Mh are the invariant masses of the two hemispheres. Here we define the 
normalized quantity 

Ml 

P^T^^ (6) 

vis 

where E^s is the total visible energy measured in hadronic events. To first order in 
perturbative QCD, and for massless partons, the heavy jet mass and thrust are related 
by r = p 
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Jet broadening measures have been proposed in Ref. In each hemisphere a, b: 

is calculated. The total jet broadening and wide jet broadening are defined by 

BT = Ba + Bh and = max(5„, 5;,), (8) 

respectively. Both Bt and Bw are identically zero in two-parton final states and 
are sensitive to the transverse structure of jets. To first order in perturbative QCD 
Bt = B\Y = ^O. 

B. Jet Rates 

Another useful method of classifying the structure of hadronic final states is in terms 
of jets. Jets may be reconstructed using iterative clustering algorithms |^ in which a 
measure i/ij, such as scaled invariant mass, is calculated for all pairs of particles i and 
j, and the pair with the smallest tjij is combined into a single particle. This procedure 
is repeated until all pairs have i/ij exceeding a value i/cut, and the jet multiplicity of the 
event is defined as the number of particles remaining. The n-jet rate RniVcut) is the 
fraction of events classified as n-jet, and the differential 2-jet rate is defined as pT 



^ R2{ycut) - R2{ycut- ^Vcut) 
D2{ycut) = T • (9) 

^Vcut 

In contrast to each event contributes to D2 at only one i/cut value. 

Several schemes have been proposed comprising different yij definitions and recom- 
bination procedures. We have applied the E, EO, P, and PO variations of the JADE 
algorithm as well as the Durham (D) and Geneva (G) schemes Q. The six defini- 



tions of the jet resolution parameter yij and recombination procedure are given below. 

In the E-scheme yij is defined as the square of the invariant mass of the pair of 
particles i and j scaled by the visible energy in the event, 

V, = (10) 

vis 
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with the recombination performed as 

Pk=Pi+Pj, (11) 

where Pi and pj are the four-momenta of the particles, and pion masses are assumed 
in calculating particle energies. Energy and momentum are explicitly conserved in this 
scheme. 

The E0-, P-, and PO-schemes are variations of the E-scheme. In the EO-scheme i/ij 
is defined by Eq. (|l3), while the recombination is defined by 

Ek = Ei + Ej (12) 

Pk = I , {Pi+Pj), (13) 

I Pi + Pj I 

where Ei and Ej are the energies and Pi and pj are the three-momenta of the particles. 
The three-momentum is rescaled so that particle k has zero invariant mass. This 
scheme does not conserve the total momentum sum of an event. 

In the P-scheme yij is defined by Eq. (0) and the recombination is defined by 

Pk = Pi + Pj (14) 
Ek = \Pk\ ■ (15) 

This scheme conserves the total momentum of an event, but does not conserve the 
total energy. 

The PO-scheme is similar to the P-scheme, but the total energy E^s in Eq- ( PI| ) is 
recalculated at each iteration according to 

Evis = X! ^k- (16) 

k 

In the D-scheme 

2rmn{Ef,E]){l-cos9,,) 
Vij = ^2 ' (17) 

^vis 

where 6ij is the angle between the pair of particles i and j. The recombination is 
defined by Eq. (pJ]). With the D-scheme a soft particle will be combined with another 
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soft particle, instead of being combined with a high-energy particle, only if the angle 
it makes with the other soft particle is smaller than the angle that it makes with the 
high-energy particle. 

The definition of Uij for the G-scheme is 

8E.E,(l-cos%) 



and the recombination is defined by Eq. (|TlD . In this scheme soft particles are combined 
as in the D-scheme. In addition, i/ij depends only on the energy of the particles to be 
combined, and not on the Ems of the event. 

C. Particle correlations 

Hadronic event observables can also be classified in terms of inclusive two-particle cor- 
relations. The energy-energy correlation {EEC) is the normalized energy-weighted 
cross section defined in terms of the angle Xij between two particles i and j in an event: 

^^^<'^' - iTTA^ S E ^Hx' - X,W, (19) 

^''events'-^X events 
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where x is an opening angle to be studied for the correlations; Ax is the angular bin 
width; and Ei and Ej are the energies of particles i and j. The angle x is taken from 
X = 0° to 180°. The shape of the EEC in the central region, x ~ 90°, is determined 
by hard gluon emission. Hadronization contributions are expected to be large in the 
collinear and back-to-back regions, x ~ 0° and 180°, respectively. The asymmetry of 
the EEC (AEEC) is defined as AEEC{x) = EEC{180° - x) - EEC{x). 

D. Angular Energy Flow 

Another procedure, related to the angle of particle emission, is to integrate the energy 
within a conical shell of opening angle x about the thrust axis. The Jet Cone Energy 
Fraction {JCEF) is defined 



as: 



■'^''''^^'> = S rj E ^Hx' - x<)dx'. (20) 

events '-^ A events i 
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where 

fpi-UT 

Xi = arccos ^ 
\ \Pi\ 

is the opening angle between a particle and the thrust axis vector, ut, whose direction 
is defined to point from the heavy jet mass hemisphere to the light jet mass hemisphere, 
and 0° < X < 180°. Hard gluon emissions contribute to the region corresponding to 
the heavy jet mass hemisphere, 90° < x ^ 180°. 




4. QCD Predictions 

The QCD predictions up to 0{a^) for all observables defined in Section 3 have the 
general form 

^ ■ ^ = ^(J/)"^ + [^(?/) + My)27rbo In/] al, (22) 

where y is the observable in question; at is the total hadronic cross section; = as/27r; 
/ = fi'^/s; Bq = (33 — 2nj)/(127r); and n/ is the number of active quark fiavors; Uf = 5 
at ^/s = Mz- We have computed the coefficients A{y) and B{y) using the EVENT 
program, which was developed by Kunszt and Nason 0]. It should be noted that a 
dependence on the QCD renormalization scale enters explicitly in the second order 
term in Eq. (P^. 

It has been found recently [jlO|, |lT], |12|, |13|, |l^, [15| that several observables, namely r. 



p, Bt, Bw, -D2(-D-scheme), and EEC, can be resummed, that is, leading and next-to- 
leading logarithmic terms can be calculated to all orders in using an exponentiation 
technique. This procedure is expected a priori to yield formulae which are less depen- 
dent on the renormalization scale. Using L = ln{l/y), the fraction R[y,as) can then 
be written in the general form 

1 fv da 

R{y,a,) = - —dy = C{as)exp{J:{as,L)} + E{y,as), (23) 
at Jo dy 
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where 

C(a,) = 1 + ^C„«:, (24) 

n=l 
oo n+1 

E(a„L) = ^a^^^G^^L-, (25) 

n=l m=l 

oo 

F{y,as) = Y.Fn{y)a:- (26) 

n=l 

The factor S to be exponentiated can be written 

E(a„ L) = L- hdasL) + fNLLia^L) + O (\ ■ (a.L)") , (27) 



oo 



where fLL^CisL) and fNLLictsL) are the leading and next-to-leading logarithms. The 



functions /ll and /nll depend only on the product a^L and are given in Refs. [|T0], [Tl 
pr2| , p!3| , |1^, |T^. The resummed calculations are thus given by an approximate expression 
for R{y, as) in the form 

R^'^-^^iy, a,) = (1 + Cia, + C^a^) exp{S™(«„ L)}, (28) 

where 

E™(a„ L) = L ■ /iL(a,L) + hiLiasL). (29) 

Whereas the leading logarithmic (L ■ /^l) and next-to-leading logarithmic {/nll) 
terms in S have been calculated, the subleading terms in Eq. (^) have not been com- 
pletely computed. However, some subleading terms included in S (Eq. (25)), as well as 
C and F, are included in the Ola^) calculation. In order to make reliable predictions, 
including hard gluon emission, with the resummed calculations it is necessary to com- 
bine them with the second order calculations, taking overlapping terms into account. 
This procedure is called matching, and four matching schemes have been proposed in 
the literature. 

The 0{al) QCD formula (Eq. (p2|)) can also be cast into the integrated form 

i?^(°^)(2/,a,) = 1+Aiy)as + Biy)al (30) 
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where A{y) and B{y) are the cumulative forms of A{y) and B(y) in Eq. p^). Taking 
the logarithm of the resummed formula (Eq. (|^)) and the 0{a1) formula (Eq. (|30D ), 



and 



\nR~{y,as) = S™(a„ L) + Cia, + - ^] + 0(aD, (31) 



ln/2^("^)(2/,a,) = A{y)a, + (B{y) - ^^^1 + 0{al). (32) 

Adding Eq. (|3T|) and Eq. (|32D , and subtracting the overlapping first and second order 
terms from Eq. (|31|), yields [|10| , p 



lni?^"+^("»)(2/,a,) = S"™(a„L) - S''™(i)(a„L) - E''"^^\a„ L) 

+Aiy)as + (^Biy) - al (33) 

where 

S™«(a„L) = Giaa.L^ + Gna.L (34) 
= G^^&lL^ + G22alL\ (35) 

Finally, one can derive j^^'^^^'^+'^i'^'i) (^y ^ q.^) by taking the exponential of Eq. (^31). This 
procedure is called InR-matching. 

In an alternative approach, the overlapping terms S'''^*"™*^^^(as, L) and S'^*^*"™'^^)(as, L) 
are subtracted from E'''^*"™(as, -^^) in the form of an exponential. The exact formula up 
to 0{af) is then obtained as follows 



i?™+'^(°'n?/,«s) = (l + Ci«, + C2«,')[exp{E™(a„L)} 

- exp{s''™W(a„L) + S'■™(2)^^^^^)y 
+l+^(2/)a,+S(?/)a^ 

= (1 + + C2«,') exp {S™(«„ L)} 

- {Cias + S™(i)(a„ L)) - [C2&I + Cia,S™(i)(a„ L) 
+ ^ {E'^™(i)(a„ L)}^ + ^) 

+ ^(2/)a, + i3(2/)a2. (36) 
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This is called R-matching, and differs from InR-matching in that the subleading term 
G2idilL is not exponentiated. In order to raise this procedure to the same level as 
the InR-matching scheme, Eq. ( ^61) may be modified by replacing S''^*"'"(as, L) and 
TJ-™^'^\as,L) with and = Gaaa^i^^ + G22tt?^^ + Gaia^L, respec- 

tively. This procedure is called modified R-matchin^ 0]. 

The predictions of these matching schemes have some troublesome features near 
the upper kinematic limit Umax because terms of third and higher order generated by 
the resummed calculations do not vanish at this limit. This situation can be corrected 
by invoking a replacement of L = y) in Eq. (^) with L' = ln{l/y — l/Vmax + !)• 
This procedure is called modified InR-matching |3^. We took the value of Vmax to be 
0.5 for T, 0.42 for p, 0.41 for Bt, 0.325 for Bw, and 0.33 for D2(D). 

Finally, in order to account for the renormalization scale dependence, fNLL{<^sL) 
should be modified to fNiLictgL) + (ctgL)^ ^If^^^T^^ f ■> ^^"^ ^iv) ^^'^ ^22 should be 
modified to B{y) + ^(?/)27r6o In / and G22 + G'i227r6o In /, respectively 0, |T5|| . 



5. Measurement of ^^(M^ 
A. Data Analysis 

The fifteen observables defined in Section 3 were calculated from the experimental 
data using charged tracks in hadronic events selected according to the criteria defined 
in Section 2. The experimental distributions D'gff){y) were then corrected for the 
effects of selection cuts, detector acceptance, efficiency, and resolution, for neutral 
particles, particle decays and interactions within the detector, and for initial state 
photon radiation, using bin-by-bin correction factors Coiy)'- 

nMC ( \, 

CD{y)r = jjMC(,A ^ (37) 



* It has also been called R-G2i-Tnatching [pO[, or intermediate matching pT[ 
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where y is the observable; i is the bin index; Dg^^j^{y)i is the content of bin i of the 
distribution obtained from reconstructed charged particles in Monte Carlo events after 
simulation of the detector; and Dhadroniy)i is that from all generated particles with 
lifetimes greater than 3 x 10^^° s in Monte Carlo events with no SLD simulation and 
no initial state radiation. The bin widths were chosen from the estimated experimental 
resolution so as to minimize bin-to-bin migration effects. The Coiy) were calculated 
using events generated with JETSET 6.3 |21| using parameter values tuned to hadronic 



e"'"e~ annihilation data In addition, the multiplicity and momentum spectra of B 
hadron decay products were tuned to T45 data [^]. The hadron level distributions are 
then given by 

DtXoAy)^ = Cnivh ■ Di%{y),. (38) 

Systematic effects were investigated using a variety of techniques. The experimental 
systematic errors arising from uncertainties in modeling the detector were estimated 
by varying the charged track and event selection criteria over wide ranges, and by 
varying the tracking efficiency and resolution in the detector simulation. In each case 
the correction factors Coiy), and hence the corrected data distributions Dhadroniy)^ 
were rederived. The data correction procedure was repeated by recalculating the cor- 
rection factors C£){y) using events generated with HERWIG 5.5 [|^. In addition, a 



matrix correction procedure was employed, in which migrations between all pairs 
of bins are accounted for individually. The differences between the data distributions 
corrected by the bin-by-bin and matrix methods were found to be much smaller than 
the statistical errors. 

The hadron level data are shown in Figs. 1-15 and listed in Tables I- VII, together 
with statistical and systematic errors; they may be compared with data from other 
experiments that have applied corrections for detector effects. The central values rep- 
resent the data corrected by the central values of the correction factors Cniy)-, which 
are shown in Figs. 16(c)-30(c). For the EEC, AEEC, and JCEF, where there are 
bin-to-bin correlations and multiple entries per event per bin, the statistical error in 
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each bin was estimated by taking the rms deviation of the contents of that bin over 50 
Monte Carlo samples, each comprising the same number of events as the data sample. 
The systematic errors derive from the uncertainties on the correction factors shown in 
Figs. 16(c)-30(c). Also shown in Figs. 1-15 are the predictions of the JETSET 7.3 
36| and HERWIG 5.5 |^ QCD + fragmentation event generators. Good agreement 



between the data and model predictions is apparent in all cases. 

Before they can be compared with the QCD predictions, the data must be corrected 
for the effects of hadronization. The correction procedure is similar to that described 
above for the detector effects. Bin-by-bin correction factors 



nMC ( \ 

■Spartan yH J 

DMC 777) 

hadron V f / 



Cnivh = T^iF^^ (39) 



where -Dpf^o„(2/)i is the content of bin i of the distribution obtained from Monte Carlo 
events generated at the parton level, were calculated and applied to the hadron level 
data distributions D'l^^^^^{y)i to obtain the parton level corrected data: 

The phenomenological hadronization models implemented in JETSET 7.3 and HER- 
WIG 5.5 were used to calculate the Cniy)- In the case of JETSET the Cniu) were 
also recalculated for values of the parton virtuality cutoff Qo |2l|, ^ in the range 



0.5 to 2.0 GeV, and for reasonable variations of the parameters All, a, and aq. The 
correction factors Cniv) are shown in Figs. 16(b)-30(b), where the bands show the 
uncertainties due to model differences and parameter variations. The parton level data 
are shown in Figs. 16(a)-30(a). The data points correspond to the central values of the 
hadronization correction factors, and the errors shown are statistical and experimental 
systematic only; the hadronization uncertainty will be considered in the next sections 
which describe the fits to determine as{M'^). 



16 



B. Measurement of as(M|) using 0{a^) calculations 

We first determined as{M'^) by comparing the Ola"^) QCD calculations for each ob- 
servable y with the corrected data at the parton level. Each calculation was fitted 
to the measured distribution 1/at ■ da/dy hj minimizing with respect to variation 
of Ajg^. In each y bin was defined using the sum in quadrature of the statistical 
and systematic errors. Fits were performed at selected values of the scale / and were 
restricted to the range in y for which the 0{a1) calculation provides a good description 
of the corrected data. 

The fit ranges in y were chosen to ensure that the parton level data and the QCD 
calculations could be compared meaningfully. The range for each observable was de- 
termined according to the following requirements: (1) the hadronization correction 
factors Cniy) satisfied 0.6 < Cniy) < 1-4; (2) the systematic uncertainties on the de- 
tector and hadronization correction factors, ACoiy) ^Cuiy) respectively, satisfied 
I ACoiy), ACniy) K 0-3; (3) three massless partons can contribute to the distribu- 
tion at 0{as) in perturbative QCD; (4) the psr degree of freedom, Xdof^ ^ at 
/ = 1 is 5.0 or less. Requirements (1) and (2) ensure that the corrected data are well 
measured and that the hadronization corrections are modeled reliably. Requirement 
(3) ensures that the kinematic regions dominated by 4-parton production at 0{a'^) are 
excluded, as the calculation is effectively leading order, and hence unreliable, in these 
regions. Requirement (4) is an empirical constraint that ensures that the QCD calcula- 
tion fits the data reasonably well; this is most relevant to exclude the so-called 'two-jet 
region' where multiple emissions of soft or collinear gluons are important and are not 
included in the O(tts) calculations, a matter discussed further in Section 5.C. These fit 
ranges are listed in Table VIII and are shown in Figs. 16-30. For illustration, fits to 
the distributions are shown in Figs. 16(a)-30(a) for the case / = 1. The data are well 
described by 0{a1) QCD within the fit ranges. Fits were also performed in the same 
ranges for different choices of the renormalization scale / such that 10~^ < / < 10^. In 
each case the fitted value of was translated to a;s(M|). The value of as(M|) 
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and the corresponding Xdof fit are shown as a function of the choice of / in 

Figs. 31-33 for all observables. 

Several features are common to the results from each observable: asiM"^) depends 
strongly on /; the fit quality is good over a wide range of /, typically / ^ 10~^, and 
there is no strong preference for a particular scale for most of the observables; at low / 
the fit quality deteriorates rapidly, and neither as(M|) nor its error can be interpreted 
meaningfully. Similar features were reported in our earlier a<((M|) measurements from 
jet rates [^] and energy-energy correlations For the oblateness the good fit region 
is / ^ 10~^, which is much higher than for the other observables. For D2 calculated in 
the E-scheme the lowest Xdof found in the region around / ~ 10~^, which is much 
lower than for the other observables. 

Figures 31-33 form a complete representation of the results of the fits of Ola"^) 
QCD to our data. It is useful, however, to quote a single value of q;s(M|), together 
with its associated uncertainties, determined from each observable. For this purpose we 
adopt the following procedure, similar to that adopted in our previous measurements 

For each observable an /-range was defined such that Xdof < ^-0 ^^'^ / — ^.0. 
The former requirement excludes the low / regions where the fit quality is poor, which 
has been shown |3^ to be due to poor convergence of the O^a^) calculations. The 
latter requirement corresponds to a reasonable physical limit /i < 2y/s. This range is 
arbitrary, but does ensure that the smallest as{M^) point (see Figs. 31(a)-33(a)) is 
considered for all variables except B^- The extrema of as(M|) values in this /-range 
were taken to define a symmetric renormalization scale uncertainty about their average, 
which we defined as the central value. The /-range, central as(M|) value, and scale 
uncertainty are listed in Table VIII for each observable. 

For most observables the statistical error on as(M|) was defined by the change in 
as(M|) corresponding to an increase in x^ of 1.0 above the lowest value within the 
/-range defined above (see Figs. 31(b)-33(b)). However, for the EEC, AEEC, and 
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JCEF, where there are strong bin-to-bin correlations, the statistical error on q;s(M|) 
was estimated by applying the same fitting procedure to ten sets of Monte Carlo events, 
each comprising the same number of events as the data sample, and taking the rms 
deviation over the ten samples. The statistical error is less than 1% of as{M^) for each 
observable, and is hsted in Table VIII. 

For each observable the experimental systematic error on q;s(M|) was estimated by 
changing the detector correction factor Cd within the systematic limits shown in Figs. 
16(c)-30(c), and by repeating the correction and fitting procedures to obtain and 
hence as{M'^) values. The systematic error, calculated from the resulting spread in 
as{Mz) values, was found to be 1-3% of ^^(Ml) for each observable and is hsted in 
Table VIII. 

For each observable the hadronization uncertainty on as(M|) was estimated by 
changing the hadronization correction factor Ch within the systematic limits shown 
in Figs. 16(b)-30(b), and by repeating the correction and fitting procedures to obtain 
Ajg^ and hence q;s(M|) values. The hadronization uncertainty, calculated from the re- 
sulting spread in q;s(M|) values, was found to be 0.4-6% of q;s(M|) for each observable 
and is hsted in Table VIII. 

The central values of as(M|) and the errors are summarized in Table IX. For each 
observable the total experimental error is the sum in quadrature of the statistical and 
experimental systematic errors, and the total theoretical uncertainty is the sum in 
quadrature of the hadronization and scale uncertainties. In all cases the theoretical 
uncertainty, which derives mainly from the scale ambiguity, dominates. This uncer- 
tainty, which arises from uncalculated higher order terms in perturbation theory, varies 
from about 3% of q;,(M|) for the AEEC to about 17% of a,(M|) for Bt- The q;s(M|) 
values from the fifteen observables are consistent within these theoretical uncertainties. 
Since the same data were used to measure all observables, and the observables are all 
highly correlated, we combine these results using an unweighted average to obtain 

q;,(M|) = 0.1225 ± 0.0026(exp.) ± 0.0109(theor.), 
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where the experimental error is the sum in quadrature of the average statistical (±0.0009) 
and average experimental systematic (±0.0024) errors, corresponding to the assump- 
tion that all are completely correlated. The theoretical error is the sum in quadrature 
of the average hadronization (±0.0024) and average scale (±0.0106) uncertainties. 

As a cross-check we combined the results by using weighted averages. Weighting 
by experimental errors yields an average as(M|) value different from the above by 
±0.0009; weighting by the total errors yields an a;s(M|) value different by —0.0013. 
These differences are of the same order as the statistical error on a single as(M|) 
measurement and are hence negligible. 

C. Measurement of as(M|) using resummed+ 0{a^) calcula- 
tions 

We next determined as{M^) by comparing the resummed±(9(a;3) calculations with the 
corrected data at the parton level for those observables for which the resummed±(9(ag) 
calculations exist, i.e. thrust (r), heavy jet mass (p), total (-Bt) and wide (Biy) 
jet broadening measures, differential 2-jet rate (-D2) calculated in the D-scheme, and 
energy-energy correlations (EEC) . We considered all four matching schemes discussed 
in Section 4, namely, InR-, modified InR-, R-, and modified R-matching. However, 
modified R-matching is not applicable to D2 because the subleading term G'2i[] is not 
calculated in this case. For the EEC InR-matching and modified InR-matching schemes 
cannot be applied reliably and were not used. 

The fit ranges were initially chosen to be the same as for the 0{a1) fits except for 
the EEC, for which the fits were performed within the angular range 90° < x ^ 154.8°, 
where the lower limit is the kinematic limit for the resummed±(9(ag) calculation. For 
the fit to D2 (D-scheme) we adopted a procedure using the matched calculation for 
0.03 < Ucut < 0.05 and the 0{a1) calculation for 0.05 < ycut < 0.33. Fits to determine 

tThe value of G21 cannot be estimated until a complete calculation of G22 is available [§9). 
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Aj^, and hence ^^(Ml), were performed as described in the previous section. For illus- 
tration Figs. 16(a)-19(a), 26(a), and 28(a) show the results of the resummed-|-(!?(Q!^) 
QCD fits using the modified InR-matching scheme with the renormalization scale factor 
/ = 1. The data are well described by the QCD calculations within the fit ranges, and 
also beyond the fit ranges into the so-called 'two-jet region' or 'Sudakov region' where 
the resummed contributions are large ||10|, |13[. This is discussed further at the end of 
this section. Figures 34-37 show (a) as(M|) and (b) the corresponding xL/; derived 
from fits at different values of /, for the four matching schemes. 

Several features should be noted from Figs. 34-37. For each matching scheme and 
each observable the dependence of as(M|) on / (Figs. 34(a)-37(a)) is weaker than 
that from the 0{al) fits (Figs. 31(a)-33(a)); the range of / for which the fit quality 
is good (Figs. 34(b)-37(b)) is in all cases smaller than the corresponding range from 
the 0{a1) fits (Figs. 31(b)-33(b)), and some observables, most notably Bt and Bw, 
do display preferences for particular scales, typically in the range 10^^ < / < 10. 
However, using the R-matching scheme we found the fit qualities for B^ and B\y to be 
very poor for all scales. For a given observable, at any given / the values of as(M|) and 
Xdo/ cire typically similar for both of the InR-matching schemes; however, the results 
from the two R-matching schemes are typically systematically different both between 
the two schemes and with respect to the two InR-matching schemes. Since there is a 
priori no strong reason to reject individual matching schemes from consideration, it is 
necessary to consider an additional theoretical uncertainty deriving from the matching 
ambiguity; this will be discussed below. 

In order to quote a single a;s(M|) value, and corresponding errors, for each ob- 
servable we applied the same procedure as for the O(a^) fits to the results from each 
matching scheme. Table X summarizes the /-ranges, central values of as(M|), and 
scale uncertainties. The experimental and hadronization systematic uncertainties were 
estimated by the methods described in Section 5.B and found to be similar to those 
from the 0{a1) analysis. For each observable we then took the average ^^(Ml) value 
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over all four matching schemes. The maximum deviation of as(M|) from the central 
value was defined as the matching uncertainty, and was added in quadrature with the 
hadronization and scale uncertainties to obtain a total theoretical uncertainty for each 
observable. The scale and matching uncertainties both derive from uncalculated higher 
order perturbative contributions and are therefore correlated, although to an unknown 
degree. The inclusion of both contributions in the total theoretical uncertainty there- 
fore represents a conservative, though not unreasonable, estimate of the effects of the 
higher order contributions. The central as(M|) value, total experimental error, defined 
as the sum in quadrature of the statistical and experimental systematic errors, and the 
total theoretical uncertainty are listed in Table XI. 

Comparing the results in Tables IX and XI it is apparent that the values of as(M|) 
from the resummed-|-(9(Q;^) fits are lower than those from the 0{al) fits by about 
3% (r), 6% (p), and 7% (Bt and Bw), but higher by about 4% (£'2(0)) and 5% 
(EEC). In addition, for all observables except D2(D), the theoretical uncertainty is 
considerably smaller for the resummed-|-(9(Q;^) case than for the 0{a1) case, despite 
the extra matching uncertainty contribution to the former. For -D2(D) the theoretical 
uncertainty is essentially the same for both 0(a;^) and resummed+(9(ag) cases, which 
may relate to the fact that the resummation of next-to-leading logarithms of tjcut to all 
orders of as is not complete^] [|T^. In all cases, however, the theoretical uncertainty is 
larger than the experimental error. 

Combining the Tesummed+0{a^) results from all six observables using an un- 
weighted average we obtain q;s(M|) = 0.1192±0.0025(exp.)±0.0070(theor.), where the 
total experimental error is the sum in quadrature of the average statistical (±0.0007) 
and average experimental systematic (±0.0024) errors, and the total theoretical error 
is the sum in quadrature of the average hadronization (±0.00 16) and average scale and 
matching(±0.0065) uncertainties. As a cross-check we combined the results by using 
weighted averages. Weighting by experimental errors yields an average as(M|) value 

■I- A complete analytic expression has recently been obtained EQ]. 
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different from the above by —0.0011; weighting by the total errors yields an q;s(M|) 
value different by —0.0015. These differences are of the same order as the statistical 
error on a single as(M|) measurement and are hence negligible. 

It is interesting to compare the resummed+0(Q;^) result with the 0{a.l) result. 
The final value quoted in Section 5.B is the average of the 0{al) results over all 15 
observables, whereas the value quoted above is the average of the resummed+0(Q;5) 
results over a subset of 6 observables. For the purposes of comparison we averaged 
the Oial) results for r, p, Bt, Bw, £'2(0), and EEC to obtain as{M^) = 0.1242 ± 
0.0026(exp.) ± 0.0132(theor.). For the same set of six observables, therefore, we find 
that the central as{M^) values derived from 0{al) and resummed+0(Q;^) fits in the 
same range of each observable are in agreement to within the (correlated) experimental 
errors, and that the theoretical uncertainty is significantly smaller when the resummed 
calculations are employed. 

From Figs. 16(a)-19(a), 26(a) and 28(a), it is clear that the resummed+0(a^) 
calculations are more successful than the 0{al) calculations in describing the two-jet 
(Sudakov) region. This implies that multiple emissions of soft gluons, which are taken 
into account in the resummed terms, contribute significantly to this region. Therefore, 
for each observable we extended the fit range into the two-jet region and extracted 
q;s(M|) as a function of the renormalization scale factor /. Requirements (l)-(3) 
(Section 5.B) were applied. In addition, for 1)2(0) we required the 5-jet production 
rate to be less than 1%; for the EEC the upper limit of the fit range was extended 
to X = 162° by applying the empirical criterion Xdof < ^- The fit ranges are listed in 
Table XII. 

The same procedure as above was applied to define a range of renormalization 
scale factor / over which to calculate a central ^^(Ml) value and scale uncertainty for 
each observable; the /-range, central as{M^) value, and scale uncertainty are hsted 
in Table XII separately for fits using each of the four matching schemes. Good fits 
with Xdof ^ could not be obtained using the R-matching scheme for r, Bt, B\y, 
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and .02(0) for any extension of the fit range beyond that used for the 0{a'^) fits. 
By comparing Tables X and XII it can be seen that the maximum change in q;s(M|) 
when the fit range is extended into the two-jet region is —0.0026 for r {InR-matching) , 
—0.0038 for p {R-matching) , —0.0009 for Bt {modified InR-matching), —0.0006 for B]y 
{modified InR-matching), —0.0045 for D2(D) {modified InR-matching), and —0.0006 
for the EEC {R-matching). These shifts are smaller than, or comparable with, the 
experimental errors, and are much smaller than the theoretical uncertainties. 

For each observable the average as{M^) value over all four matching schemes, and 
the matching uncertainty, were calculated as before. The central q;s(M|) value, the 
total experimental error, and the total theoretical uncertainty, defined as before, are 
listed in Table XIII. Averaging over the six observables, as above, then yields 

q;,(M|) = 0.1181 ± 0.0024(exp.) ± 0.0057(theor.), 

which is in good agreement with the above average of results from the restricted fit 
ranges. 

6. Conclusions 

We have measured the strong coupling q;s(M|) by analyses of fifteen different observ- 
ables that describe the hadronic final states of about 60,000 decays recorded by the 
SLD experiment. The observables comprise six event shapes (r, p, Bt, B^, O, and 
C), differential 2-jet rates (-D2) defined by six different jet resolution/recombination 
schemes (E, EO, P, PO, D, and G), energy-energy correlations {EEC) and their asym- 
metry {AEEC), and the jet cone energy fraction {JCEF). The quantity JCEF has 
been measured for the first time. Our measured distributions of these observables are 
reproduced by the JETSET and HERWIG Monte Carlo simulations of hadronic Z° 
decays. The coupling was determined by fitting perturbative QCD calculations to the 
data corrected to the parton level. Perturbative QCD calculations complete to 0{a1) 
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were used for all 15 observables. In addition, recently-performed resummed calcula- 
tions were matched to the 0{a1) calculations using four matching schemes and apphed 
to the six observables for which the resummed calculations are available. 

We find that the 0{a1) calculations are able to describe the data in the hard 3-jet 
region of all 15 observables for a wide range of the QCD renormalization scale factor 
/. The fitted as(M|) value depends strongly both on the choice of /, which limits 
the precision of the as(M|) measurement from each observable, and on the choice 
of observable. The AEEC shows the smallest renormalization scale uncertainty of 
about 3%, which is just larger than the experimental error. The ^^(Mf ) values deter- 
mined from jet rates and energy-energy correlations are consistent with our previous 
measurements Q within experimental errors. The as(M|) values from the various 
observables are consistent with each other only within the scale uncertainties. The large 
scale uncertainties and systematically different q;s(M|) values determined from differ- 
ent observables imply that the uncalculated 0{al) perturbative QCD contributions 
are significant and cannot be ignored if as(M|) is to be determined with a precision 
of better than 10%. 

The resummed+(9(a^) calculations yield a reduced renormalization scale depen- 
dence of q;s(M|), and fit a wider kinematic region, including the two-jet or Sudakov 
region, and give similar fitted values of q;s(M|) to the 0{a1) case. However, the differ- 
ent matching schemes give different q;s(M|) values, which reflects a residual uncertainty 
in the inclusion of terms in the resummed-|-(9(Q;^) calculations. For all observables ex- 
cept -D2(D) the theoretical uncertainty is smaller than in the 0{al) case, but still 
dominates the uncertainty in the measurement of as(M|). Again, the as(M|) values 
derived from jet rates and energy-energy correlations are consistent with our previous 
measurements p, ^ within experimental errors, and the values determined from the 
six observables are consistent within theoretical uncertainties. 

Figure 38 summarizes the measured q;s(M|) values from all fifteen observables using 
calculations, and from the six observables using Tesummed+0{af) calculations 
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in the extended kinematic region. Since the same data were used to measure all ob- 
servables, and the observables are highly correlated, we combined the results by taking 
unweighted averages of the q;s(M|) values and experimental and theoretical errors, 
obtaining 

a^(M|) = 0.1225 ± 0.0026(exp.) ±0.0109(theor.) 0{a^J 

as{Ml) = 0.1181 ±0.0024(exp.) ±0.0057(theor.) resummed + C(a^), 

where in both cases the theoretical uncertainty is dominated by the lack of knowledge of 
higher order terms in the QCD calculations. Our estimate of the theoretical uncertainty 
is larger than that quoted by some of the LEP experiments because we have considered 
more observables and wider variations of the renormalization scale, and have taken 
unweighted averages. These average values are shown in Fig. 38; they are consistent 
with measurements from other e"'"e~ experiments at the resonance ^ |3^, ^ 



and from lower energy e e~ and deep inelastic scattering experiments ||42| . 

One expects a priori the as(M|) value determined from a resummed+(9(a^) fit to 
be more reliable than that from an 0(q;^) fit. However, the former is only available for 
six of the fifteen observables. In order to quote a final result, therefore, we took the 
unweighted average of the as(M|) values and uncertainties over the combined set of six 
resummed+(9(a^) results and nine 0{a1) results for which there is no corresponding 
resummed+C(a^) result. This yields a final average of 

a,(M|) = 0.1200 ± 0.0025(exp.) ± 0.0078(theor.), 

also shown in Fig. 38, corresponding to = 253l^|6 MeV. 
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T 


j^^±{stat.)±iexp. sys.) 


P 


^^±istat.)±{exp. sys.) 





- 


0.02 


7.01 ±0.10 


± 0.50 


0.0 


- 0.02 


10.53 ±0.12 ±0.41 





02- 





04 


16.10 ±0.15 


±0.15 


0.02 


- 0.04 


17.38 ±0.15 ±0.14 





04- 





06 


8.67 ±0.11 


±0.05 


0.04 


- 0.08 


6.21 ±0.07 ±0.16 





06- 





08 


5.08 ±0.08 


±0.16 


0.08 


- 0.12 


2.39 ±0.04 ±0.09 





08- 





12 


2.91 ±0.04 


±0.06 


0.12 


- 0.18 


1.08 ±0.02 ±0.04 





12- 





16 


1.57 ±0.03 


±0.05 


0.18 


- 0.24 


0.404 ±0.014 ±0.021 





16- 





20 


0.917 ±0.025 


± 0.028 


0.24 


- 0.32 


0.102 ±0.006 ±0.010 





20- 





26 


0.495 ±0.015 


± 0.025 


0.32 


- 0.40 


0.0047 ±0.0013 ±0.0008 





26- 





32 


0.227 ±0.010 


±0.016 











32- 





38 


0.061 ±0.005 


±0.006 











38- 





44 


0.003 ±0.001 


± 0.003 









Table I. Distributions of r and p (see text). The data were corrected for detector 
effects and for initial state photon radiation. The first error is statistical, and the 
second represents the experimental systematic uncertainty. 
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i;m^±istat.)±{exp. sys.) 


Bw 


± ±(earp. sys.) 


0.0 


-0.02 


0.018 ± 0.005 ±0.007 


0.0 


-0.02 


0.570 ±0.028 ±0.213 


0.02 


- 0.04 


1.36 ±0.04 ±0.18 


0.02 


- 0.04 


13.86 ±0.14 ±0.45 


0.04 


- 0.06 


8.81 ±0.11 ±0.32 


0.04 


-0.06 


11.71 ±0.13 ±0.20 


0.06 


-0.08 


10.64 ±0.12 ±0.16 


0.06 


-0.08 


7.38 ±0.10 ±0.11 


0.08 


-0.12 


6.52 ±0.07 ±0.10 


0.08 


-0.12 


4.29 ±0.05 ±0.08 


0.12 


- 0.16 


3.65 ±0.05 ±0.04 


0.12 


- 0.16 


2.185 ±0.038 ±0.128 


0.16 


- 0.20 


2.10 ±0.04 ±0.06 


0.16 


- 0.20 


1.12 ± 0.028 ±0.061 


0.20 


- 0.26 


1.12 ±0.02 ±0.03 


0.20 


- 0.26 


0.403 ±0.014 ±0.025 


0.26 


-0.32 


0.384 ±0.013 ±0.023 


0.26 


-0.32 


0.403 ±0.014 ±0.025 


0.32 


-0.38 


0.050 ±0.005 ±0.011 


0.32 


-0.38 


0.030 ±0.004 ±0.005 



Table II. Distributions of and Bw (see text). The data were corrected for detector 
effects and for initial state photon radiation. The first error is statistical, and the 
second represents the experimental systematic uncertainty. 



o 


j^j^±{stat.)±iexp. sys.) 


C 


j;^±{stat.)±{exp. sys.) 


0.0 


-0.02 


9.07 ±0.11 ±0.19 








-0.04 


0.166 ±0.011 ±0.015 


0.02 


- 0.04 


11.28 ±0.12 ±0.20 





04 


-0 


08 


1.76 ±0.03 ±0.04 


0.04 


- 0.08 


5.98 ±0.06 ±0.07 





08 


-0 


12 


4.01 ±0.05 ±0.09 


0.08 


- 0.12 


3.16 ±0.05 ±0.06 





12 


-0 


18 


3.57 ±0.04 ±0.10 


0.12 


- 0.18 


1.77 ±0.03 ±0.03 





18 


-0 


24 


2.30 ±0.03 ±0.02 


0.18 


- 0.24 


0.935 ±0.021 ±0.028 





24 


- 


32 


1.54 ±0.02 ±0.016 


0.24 


- 0.32 


0.523 ±0.013 ±0.013 





32 


- 


40 


1.07 ±0.02 ±0.03 


0.32 


- 0.40 


0.223 ±0.009 ±0.010 





40 


- 


52 


0.718 ±0.013 ±0.024 


0.40 


- 0.50 


0.052 ±0.004 ±0.003 





52 


- 


64 


0.491 ±0.011 ±0.013 











64 


- 


76 


0.311 ±0.008 ±0.022 











76 


- 


88 


0.146 ±0.006 ±0.012 











88 


- 1 





0.012 ±0.002 ±0.001 



Table III. Distributions of O and C (see text). The data were corrected for detector 
effects and for initial state photon radiation. The first error is statistical, and the 
second represents the experimental systematic uncertainty. 
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E-scheme 


EO-scheme 


P-scheme 


Vcut 


DiiVcut) ± {stat.) ± {exp. sys.) 


D^iVcut) ± {stat.) ± {exp. sys.) 


D2{ycut) ± {stat.) ± {exp. sys.) 


0.005 


0.669 ±0.060 ±0.080 


28.95 ±0.39 ± 1.44 


41.80 ±0.47 ±2.43 


0.010 


2.60 ±0.12 ±0.12 


25.25 ±0.37 ±0.50 


31.06 ±0.41 ±0.63 


0.015 


7.07 ±0.20 ±0.27 


19.93 ±0.33 ±0.53 


21.24 ±0.34 ±0.28 


0.02 


10.48 ±0.24 ±0.66 


15.85 ±0.29 ±1.04 


14.96 ±0.28 ±0.54 


0.03 


12.28 ±0.18 ±0.39 


11.66 ±0.18 ±0.15 


10.82 ±0.17 ±0.37 


0.05 


10.89 ±0.12 ±0.34 


7.01 ±0.10 ±0.19 


6.35 ±0.09 ±0.23 


0.08 


7.22 ±0.08 ±0.22 


3.85 ±0.06 ±0.05 


3.16 ±0.05 ±0.09 


0.12 


3.81 ±0.05 ±0.11 


2.02 ±0.04 ±0.07 


1.61 ±0.03 ±0.08 


0.17 


1.97 ±0.03 ±0.05 


1.08 ±0.02 ±0.04 


0.791 ± 0.021 ± 0.037 


0.22 


0.987 ±0.023 ±0.034 


0.537 ±0.017 ±0.026 


0.317 ±0.013 ±0.024 


0.28 


0.467 ±0.015 ±0.017 


0.204 ±0.010 ±0.015 


0.069 ±0.006 ±0.005 


0.33 


0.178 ±0.009 ±0.024 


0.068 ±0.006 ±0.021 


0.008 ± 0.002 ± 0.007 



Table IV. -D2(l/c«t) calculated in the E-scheme, the EO-scheme, and the P-scheme (see 
text). The data were corrected for detector effects and for initial state photon radiation. 
The first error is statistical, and the second represents the experimental systematic 
uncertainty. 
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PO-scheme 


D-scheme 


G-scheme 


Vcut 


DiiVcut) ± {stat.) ± {exp. sys.) 


D2{ycut) ± {stat.) ± {exp. sys.) 


D2{ycut) ± {stat.) ± {exp. sys.) 


0.005 


39.78 ±0.46 ±2.41 


101.06 ±0.74 ±2.29 


7.67 ±0.20 ±1.01 


0.010 


29.85 ±0.40 ±0.78 


26.85 ±0.38 ±0.34 


33.63 ±0.43 ±0.84 


0.015 


20.49 ±0.33 ±0.36 


14.13 ±0.28 ±0.40 


31.71 ±0.41 ± 1.01 


0.02 


14.52 ±0.28 ±0.23 


9.00 ±0.22 ±0.44 


20.46 ±0.33 ±0.55 


0.03 


10.65 ±0.17 ±0.37 


6.02 ±0.13 ±0.17 


11.71 ±0.18 ±0.20 


0.05 


6.36 ±0.09 ±0.19 


3.30 ±0.07 ±0.11 


5.55 ±0.09 ±0.12 


0.08 


3.21 ±0.05 ±0.12 


1.66 ±0.04 ±0.07 


3.20 ±0.05 ±0.06 


0.12 


1.64 ±0.03 ±0.07 


0.831 ±0.024 ±0.038 


1.92 ±0.04 ±0.05 


0.17 


0.944 ±0.023 ±0.057 


0.406 ±0.015 ±0.033 


1.25 ±0.03 ±0.03 


0.22 


0.433 ±0.015 ±0.038 


0.173 ±0.010 ±0.011 


0.768 ± 0.020 ± 0.027 


0.28 


0.169 ±0.009 ±0.015 


0.084 ±0.006 ±0.013 


0.409 ±0.014 ±0.019 


0.33 


0.034 ±0.004 ±0.008 


0.027 ±0.004 ±0.048 


0.111 ±0.007 ±0.018 



Table V. D2{ycut) calculated in the PO-scheme, the D-scheme, and the G-scheme (see 
text). The data were corrected for detector effects and for initial state photon radiation. 
The first error is statistical, and the second represents the experimental systematic 
uncertainty. 
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X (deg.) 



EEC{rad-'^) ±{stat.) ±{exp. sys.) 



X (deg.) 



EEC{rad-^) ±{stat.) ±{exp. sys.) 



0.0 
3.6 
7.2- 
10.8- 
14.4- 
18.0 
21.6 
25.2- 
28.8- 
32.4- 
36.0 
39.6 
43.2- 
46.8- 
50.4- 
54.0- 
57.6- 
61.2 
64.8- 
68.4- 
72.0- 
75.6 
79.2- 
82.8- 
86.4- 



3.6 

7.2 

10.8 

14.4 

18.0 

21.6 

25.2 

■28.8 
32.4 
36.0 
39.6 
43.2 
46.8 
50.4 
54.0 
57.6 
61.2 
64.8 
68.4 
72.0 
75.6 
79.2 

■82.8 
86.4 
90.0 



2.265 
1.316 
0.874 
0.598 
0.425 
0.310 
0.241 
0.199 
0.168 
0.146 
0.128 
0.118 
0.1099 
0.1014 
0.0935 
0.0901 
0.0867 
0.0827 
0.0802 
0.0764 
0.0770 
0.0752 
0.0736 
0.0751 
0.0744 



±0.006 

±0.006 

±0.004 

±0.003 

±0.002 

±0.002 

±0.001 

±0.001 

±0.001 

±0.001 

±0.001 

±0.001 

±0.0008 

±0.0009 

±0.0008 

±0.0009 

±0.0008 

±0.0009 

±0.0010 

±0.0009 

±0.0010 

±0.0008 

±0.0008 

±0.0010 

±0.0010 



±0.055 

±0.032 

±0.020 

±0.019 

±0.011 

±0.014 

±0.005 

±0.005 

±0.006 

±0.005 

±0.004 

±0.003 

±0.0026 

±0.0031 

±0.0027 

±0.0021 

±0.0023 

±0.0023 

±0.0018 

±0.0031 

±0.0010 

±0.0031 

±0.0013 

±0.0015 

±0.0014 



90.0 
93.6 
97.2- 
100.8 
104.4 
108.0 
111.6 
115.2 
118.8- 
122.4 
126.0 
129.6 
133.2 
136.8 
140.4 
144.0 
147.2 
151.2 
154.8 
158.4 ■ 
162.0 
165.6 
169.2 
172.8 
176.4 



93.6 

97.2 

100.8 

104.4 

108.0 

111.6 

115.2 

118.8 

122.4 

126.0 

129.6 

133.2 

136.8 

140.4 

144.0 

147.2 

151.2 

154.8 

158.4 

162.0 

165.6 

169.2 

172.8 

176.4 

180.0 



0.0761 
0.0764 
0.0777 
0.0809 
0.0834 
0.0874 
0.0931 
0.0968 
0.1030 
0.111 
0.121 
0.136 
0.151 
0.170 
0.193 
0.225 
0.265 
0.320 
0.390 
0.491 
0.636 
0.847 
1.098 
1.276 
0.764 



±0.0009 

±0.0009 

±0.0009 

±0.0012 

±0.0010 

±0.0010 

±0.0013 

±0.0012 

±0.0012 

±0.001 

±0.001 

±0.002 

±0.002 

±0.002 

±0.002 

±0.002 

±0.002 

±0.003 

±0.003 

±0.003 

±0.004 

±0.006 

±0.005 

±0.007 

±0.007 



±0.0013 

±0.0025 

±0.0023 

±0.0016 

±0.0024 

±0.0022 

±0.0015 

±0.0038 

±0.0070 

±0.002 

±0.007 

±0.003 

±0.004 

±0.005 

±0.006 

±0.008 

±0.007 

±0.008 

±0.013 

±0.017 

±0.012 

±0.007 

±0.009 

±0.044 

±0.050 



Table VI. The EEC (see text). The data were corrected for detector effects and for 
initial state photon radiation. The first error is statistical, and the second represents 
the experimental systematic uncertainty. 
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X (deg.) 


AEEC{rad-^) ±{stat.) ±{exp. sys.) 


X (deg.) 


JCEF{rad-^) ±{stat.) ±{exp. sys.) 


0.0 - 


3.6 




90.0- 


93.6 


0.0274 ±0.0016 


±0.0010 


3.6 - 


7.2 




93.6- 


97.2 


0.0403 ±0.0020 


±0.0012 


7.2- 


10.8 


0.224 ±0.010 ±0.002 


97.2- 


100.8 


0.0442 ±0.0026 


±0.0010 


10.8- 


14.4 


0.249 ±0.009 ±0.005 


100.8 - 


104.4 


0.0523 ±0.0029 


±0.0023 


14.4- 


18.0 


0.211 ±0.006 ±0.005 


104.4 - 


108.0 


0.0566 ±0.0029 


±0.0024 


18.0 - 


21.6 


0.181 ±0.004 ±0.005 


108.0- 


111.6 


0.0613 ±0.0034 


±0.0026 


21.6 - 


25.2 


0.148 ±0.004 ±0.006 


111.6- 


115.2 


0.0725 ±0.0039 


±0.0017 


25.2- 


28.8 


0.121 ±0.003 ±0.004 


115.2- 


118.8 


0.0832 ±0.0055 


±0.0046 


28.8- 


32.4 


0.0972 ±0.0024 ±0.0029 


118.8- 


122.4 


0.0858 ±0.0051 


±0.0016 


32.4- 


36.0 


0.0785 ±0.0022 ±0.0062 


122.4- 


126.0 


0.0944 ±0.0043 


±0.0024 


36.0- 


39.6 


0.0645 ±0.0017 ±0.0024 


126.0- 


129.6 


0.1051 ±0.0061 


±0.0055 


39.6 - 


43.2 


0.0513 ±0.0020 ±0.0026 


129.6 - 


133.2 


0.114 ±0.005 


±0.002 


43.2- 


46.8 


0.0413 ±0.0015 ±0.0027 


133.2- 


136.8 


0.131 ±0.005 


±0.005 


46.8- 


50.4 


0.0346 ±0.0016 ±0.0021 


136.8 - 


140.4 


0.148 ±0.005 


±0.006 


50.4- 


54.0 


0.0275 ±0.0013 ±0.0060 


140.4 - 


144.0 


0.169 ±0.007 


±0.004 


54.0- 


57.6 


0.0213 ±0.0010 ±0.0024 


144.0- 


147.2 


0.188 ±0.007 


±0.005 


57.6- 


61.2 


0.0163 ±0.0008 ±0.0073 


147.2- 


151.2 


0.228 ±0.008 


±0.009 


61.2- 


64.8 


0.0141 ±0.0007 ±0.0026 


151.2- 


154.8 


0.275 ±0.009 


±0.010 


64.8- 


68.4 


0.0129 ±0.0010 ±0.0008 


154.8 - 


158.4 


0.329 ±0.011 


±0.013 


68.4- 


72.0 


0.0110 ±0.0007 ±0.0025 


158.4- 


162.0 


0.414 ±0.011 


±0.019 


72.0- 


75.6 


0.0064 ±0.0005 ±0.0017 


162.0- 


165.6 


0.551 ±0.012 


±0.013 


75.6 - 


79.2 


0.0058 ±0.0006 ±0.0029 


165.6- 


169.2 


0.751 ±0.021 


±0.021 


79.2- 


82.8 


0.0041 ±0.0004 ±0.0020 


169.2- 


172.8 


1.095 ±0.024 


±0.019 


82.8- 


86.4 


0.0012 ±0.0002 ±0.0038 


172.8 - 


176.4 


1.639 ±0.032 


±0.034 


86.4- 


90.0 


0.0017 ±0.0008 ±0.0016 


176.4- 


180.0 


1.530 ±0.039 


±0.049 



Table VII. The AEEC and JCEF (see text). The data were corrected for detector 
effects and for initial state photon radiation. The first error is statistical, and the 
second represents the experimental systematic uncertainty. 
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observable 


fit range 


/-range 


a,(M|) 


uuccrt Mint icH 


Stat. 


exp. sys. 




scale 


T 


0.06 - 0.32 


2 X 10~^ — 4 


0.1245 


±0.0008 


±0.0017 


±0.0026 


±0.0201 


P 


0.04 — 0.32 


1.5 X 10 — 4 


0.1273 


itU.UUUo 


-i-H nn9n 

itU.uUZU 


ztU.UUUO 


ztu.uuyo 




n 1 9 n 
u. iz — u.oz 


0. / X iU — 4 


n 1 979 


+0 0008 


+0 0020 


+0 0033 


+0 0220 


By/ 


Ofi — 2fi 


2 X 10~^ — 4 


0.1196 


±0.0008 


±0.0026 


±0.0024 


±0.0072 


o 


0.08-0.32 


2 X 10-1 -4 


0.1343 


±0.0013 


±0.0015 


±0.0087 


±0.0082 


c 


0.24-0.76 


4 X 10-4-4 


0.1233 


±0.0009 


±0.0019 


±0.0032 


±0.0186 


£'2(E) 


0.08-0.28 


5 X lO-'"^ - 4 


0.1273 


±0.0006 


±0.0016 


±0.0022 


±0.0217 


D2(E0) 


0.05-0.28 


1.2 X 10-2 _ 4 


0.1175 


±0.0007 


±0.0027 


±0.0010 


±0.0083 


i^2(P) 


0.05-0.22 


5.5 X 10-3 -4 


0.1207 


±0.0008 


±0.0033 


±0.0025 


±0.0053 


^2(P0) 


0.05 - 0.28 


1.2 X 10-2 -4 


0.1190 


±0.0009 


±0.0031 


±0.0020 


±0.0057 


i?2(D) 


0.03-0.22 


1.7 X 10-3 -4 


0.1245 


±0.0011 


±0.0032 


±0.0007 


±0.0077 


I>2(G) 


0.12-0.28 


4 X 10-3-4 


0.1191 


±0.0008 


±0.0014 


±0.0029 


±0.0043 




36.0° - 154.8° 


3.5 X 10-3-4 


0.1222 


±0.0008 


±0.0030 


±0.0021 


±0.0121 


AEEC 


18.0° - 68.4° 


9 X 10-2-4 


0.1121 


±0.0012 


±0.0032 


±0.0017 


±0.0031 


JCEF 


100.8° - 158.4° 


5 X 10-3 - 4 


0.1185 


±0.0007 


±0.0027 


±0.0008 


±0.0045 



Table VIII. Observables used in Oic^^ QCD fits. For each the fit range, the range 
of the renormalization scale factor considered, central o.aiM'z) value, statistical and 
experimental systematic errors, and hadronization and scale uncertainties are shown. 
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observable 




cxp. error 


theoretical uncertainty 


T 


0.1245 


±0.0019 


±0.0203 


P 


0.1273 


±0.0022 


±0.0096 


Bj" 


0.1272 


±0.0022 


±0.0222 


Bw 


0.1196 


±0.0027 


±0.0076 


O 


0.1343 


±0.0020 


±0.0120 


c 


0.1233 


±0.0021 


±0.0189 


D2(E) 


0.1273 


±0.0017 


±0.0218 


D2(E0) 


0.1175 


±0.0028 


±0.0084 


D2(P) 


0.1207 


±0.0034 


±0.0059 


£'2(P0) 


0.1190 


±0.0032 


±0.0060 


£>2(D) 


0.1245 


±0.0034 


±0.0077 


£>2(G) 


0.1191 


±0.0016 


±0.0052 




0.1222 


±0.0031 


±0.0123 




0.1121 


±0.0034 


±0.0035 


JCEF 


0.1185 


±0.0028 


±0.0046 



Table IX. The q;,(M|) values derived from 0{al) QCD fits. 
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observable 


fit range 


Ini? matching 


mod. Ini? matctiing 


R matching 


mod. R matching 


a,(M|)±Aa, 
/-range 


a,(M|)± Aa, 
/-range 


a,(M|)± Aa, 
/-range 


as{Ml)±Aas 
/-range 


T 


0.06 - 0.32 


0.1196 ±0.0089 
2.7 X 10-3-4 


0.1203 ±0.0089 
2.7 X 10-3 - 4 


0.1226 ±0.0110 
1.9 X 10-3-4 


0.1187 ±0.0091 
2.3 X 10-3 - 4 


P 


0.04 - 0.32 


0.1151 ±0.0039 
1.1 X 10-2-4 


0.1162 ±0.0047 
1.1 X 10-2-4 


0.1178 ±0.0061 
4.9 X 10-3-4 


0.1146 ±0.0044 
1.0 X 10-2-4 


Bx 


0.12-0.32 


0.1175 ±0.0030 
6.7 X 10-2-4 


0.1211 ±0.0015 
3.0 X 10-1-4 




0.1177± 0.0017 
3.6 X 10-2-4 


Bw 


0.06 - 0.26 


0.1083 ±0.0016 
8.2 X 10-2-4 


0.1095 ±0.0003 
1.9 X 10-1-4 




0.1107 ±0.0034 
4.9 X 10-2-4 


£'2(D) 


0.03 - 0.22 


0.1312 ±0.0060 
1.5 X 10-1-4 


0.1313 ±0.0059 
1.6 X 10-1-4 


0.1251 ±0.0053 
7.0 X 10-2 - 4 


N/A 


EEC 


90.0° - 154.8° 


N/A 


N/A 


0.1239 ±0.0049 
6.1 X 10-2-4 


0.1336 ±0.0028 
2.7 X 10-1-4 



Table X. Observables used in resummed+0(Q;5) fits. For each the fit range, the range 
of the renormahzation scale factor considered, the central q;s(M|) value, and scale 
uncertainty (Acts) are given. Results are shown separately for each of the four matching 
schemes considered. Acceptable fits to the data could not be obtained for Bt and Bw 
with the R-matching scheme. 
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observable 


as{Ml) 


exp. error 


theoretical uncertainty 


T 


0.1180 


±0.0018 


±0.0115 


P 


0.1163 


±0.0020 


±0.0064 




0.1160 


±0.0020 


±0.0048 


Bw 


0.1074 


±0.0025 


±0.0042 


^2(D) 


0.1297 


±0.0035 


±0.0073 


EEC 


0.1279 


±0.0032 


±0.0069 



Table XL The ^^(Ml) values derived from resummed+C(Q;^) QCD fits. 
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observable 


fit range 


Ini? matching 


mod. Ini? matcliing 


R matcliing 


mod. R matching 


as{Ml)±Aas 
/-range 


as{Ml) ± Aas 
/-range 


as(M|)± Aa^ 
/-range 


as{Ml) ± Aa^ 
/-range 


T 


0.02 - 0.32 


0.1170 ±0.0086 
7.0 X 10-2 - 4 


0.1184 ±0.0075 
1.4 X 10-1-4 




0.1191 ±0.0045 
6.3 X 10-1 -4 


P 


0.02 - 0.32 


0.1153 ±0.0071 
2.6 X 10-2-4 


0.1146 ±0.0072 
3.4 X 10-2 - 4 


0.1140 ±0.0054 
2.0 X 10-1-4 


0.1124 ±0.0071 
4.0 X 10-2 - 4 


Bt 


0.04 - 0.32 


0.1177 ±0.0040 
2.0 X 10-1-4 


0.1202 ±0.0021 
6.7 X 10-2 - 4 




0.1175 ±0.0023 
1.1 X 10-1 -4 


Bw 


0.04 - 0.26 


0.1078 ±0.0024 
1.4 X 10-1-4 


0.1089 ±0.0014 
2.8 X 10-1-4 




0.1106 ±0.0032 
5.4 X 10-2-4 


£'2(D) 


0.01-0.22 


0.1269 ± 0.0026 

1.3 X 10^1-4 


0.1268 ± 0.0025 
1.3 X 10-1-4 




N/A 


EEC 


90.0° - 162.0° 


N/A 


N/A 


0.1233 ±0.0043 
6.9 X 10-2-4 


0.1337 ± 0.0027 
5.0 X 10-1-4 



Table XII. Observables used in resummed+C(Q;5) fits with the fit ranges extended into 
the two-jet region. For each the fit range, the range of the renormalization scale factor 
considered, the central q;s(M|) value, and scale uncertainty (Aa^) are given. Results 
are shown separately for each of the four matching schemes considered. Acceptable 
fits to the data could not be obtained for r, Bt, Bw, and £^2(0) with the R-matching 
scheme. 
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observable 




exp. error 


theoretical uncertainty 


r 


0.1159 


±0.0017 


±0.0090 


P 


0.1144 


±0.0019 


±0.0074 


Bj" 


0.1157 


±0.0020 


±0.0053 


Bw 


0.1070 


±0.0025 


±0.0041 


£>2(D) 


0.1274 


±0.0034 


±0.0027 


EEC 


0.1285 


±0.0032 


±0.0068 



Table XIII. The as(M|) values derived from resummed+0(Q;^) QCD fits with the fit 
ranges extended into the two- jet region. 
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Figure captions 



Figure 1. The measured thrust distribution corrected to the hadron level. The error 

bars include the statistical and experimental systematic errors added in quadrature. 

The curves show the predictions of the QCD parton shower models JETSET 7.3 (solid 

line) and HERWIG 5.5 (dashed line). 

Figure 2. As Fig. 1 but for the heavy jet mass. 

Figure 3. As Fig. 1 but for the total jet broadening. 

Figure 4. As Fig. 1 but for the wide jet broadening. 

Figure 5. As Fig. 1 but for the oblateness. 

Figure 6. As Fig. 1 but for the C-parameter. 

Figure 7. As Fig. 1 but for the differential 2-jet rate with the E-scheme. 
Figure 8. As Fig. 1 but for the differential 2-jet rate with the EO-scheme. 
Figure 9. As Fig. 1 but for the differential 2-jet rate with the P-scheme. 
Figure 10. As Fig. 1 but for the differential 2-jct rate with the PO-schcmc. 
Figure 11. As Fig. 1 but for the differential 2-jet rate with the D-scheme. 
Figure 12. As Fig. 1 but for the differential 2-jet rate with the G-scheme. 
Figure 13. As Fig. 1 but for the energy-energy correlation EEC. 
Figure 14. As Fig. 1 but for the asymmetry of the energy-energy correlation AEEC. 
Figure 15. As Fig. 1 but for the jet cone energy fraction JCEF. 
Figure 16. (a) The measured thrust distribution corrected to the parton level. The 
error bars include the statistical and experimental systematic errors added in quadra- 
ture. The curves show the predictions of the 0(q;^) calculations (sohd hue) and the 
resummed-|-C(Q;f ) calculations with modified InR-matching (dashed line). The renor- 
malization scale factor was fixed to 1. Sizes of the (b) hadronization correction and (c) 
detector correction factors; the widths of the bands indicate the systematic uncertain- 
ties. 

Figure 17. As Fig. 16 but for the heavy jet mass. 
Figure 18. As Fig. 16 but for the total jet broadening. 
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Figure 


19. 


As 


Fig. 


16 but for the wide jet broadening. 


Figure 


20. 


As 


Fig. 


16 but for the oblateness. 


Figure 


21. 


As 


Fig. 


16 but for the C-parameter. 


Figure 


22. 


As 


Fig. 


16 but for the differential 2-jet rate with the E-scheme. 


Figure 


23. 


As 


Fig. 


16 but for the differential 2-jet rate with the EO-scheme. 


Figure 


24. 


As 


Fig. 


16 but for the differential 2-jet rate with the P-scheme. 


Figure 


25. 


As 


Fig. 


16 but for the differential 2-jet rate with the PO-scheme. 


Figure 


26. 


As 


Fig. 


16 but for the differential 2-jet rate with the D-scheme. 


Figure 


27. 


As 


Fig. 


16 but for the differential 2-jet rate with the G-scheme. 


Figure 


28. 


As 


Fig. 


16 but for the energy-energy correlation EEC. 


Figure 


29. 


As 


Fig. 


16 but for the asymmetry of the energy-energy correlation AEEC. 


Figure 


30. 


As Fig. 


16 but for the jet cone energy fraction JCEF. 


Figure 


31. 


(a 


) q;s(M|) and (b) Xdof from the 0{a'l) fits to the event shapes as a 



function of renormalization scale factor / (see text). 

Figure 32. (a) as{Mz) and (b) Xdof from the 0{af) fits to the jet rates as a function 
of renormalization scale factor / (see text). 

Figure 33. (a) ^^(Ml) and (b) Xdof from the 0{a1) fits to the particle correlations 

and angular energy flow as a function of renormalization scale factor / (sec text). 

Figure 34. (a) q;s(M|) and (b) Xdof from the resummed-|-(9(Q;g) fits with InR- 

matching as a function of renormalization scale factor / (see text) . 

Figure 35. (a) q;s(M|) and (b) Xdoj from the resummed-|-0(Q;^) fits with modified InR- 

matching as a function of renormalization scale factor / (see text). 

Figure 36. (a) as{M'^) and (b) Xdof from the resummcd+(9(a;g) fits with R-matching 

as a function of renormalization scale factor / (see text). The Xdof values for Bt and 

Bw are larger than 10 for all /. 

Figure 37. (a) as{Mz) and (b) Xdof from the resummed-|-0(Q;^) fits with modified R- 

matching as a function of renormalization scale factor / (sec text). 

Figure 38. Compilation of final values of as{M^). For each observable the solid 
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bar denotes the experimental error, while the dashed bar shows the total uncertainty 
comprising the experimental error and theoretical uncertainty in quadrature. Shown 
separately for the 0{al) results and resummed+(9(Q;^) results arc a vertical line and 
a shaded region representing the average ag^M^) value and uncertainty, respectively, 
in each case. Also shown is the final average of six resummed+(9(Q;^) and nine 0{a^g) 
results indicated by stars. 
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